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Abstract 

The Pauli operator describes the energy of a nonrelativistic quantum particle with 
spin ^ in a magnetic field and an external potential. Bounds on the sum of the negative 
eigenvalues are called magnetic Lieb-Thirring (MLT) inequalities. The purpose of this 
paper is twofold. First, we prove a new MLT inequality in a simple way. Second, we give 
a short summary of our recent proof of a more refined MLT inequality ES-Py] and we 
explain the differences between the two results and methods. The main feature of both 
estimates, compared to earlier results, is that in the large field regime they grow with 
the optimal (first) power of the strength of the magnetic field. As a byproduct of the 
method, we also obtain optimal upper bounds on the pointwise density of zero energy 
eigenfunctions of the Dirac operator. 
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1 Introduction 

1.1 Magnetic Lieb-Thirring inequalities 

Since the seminal paper of Lieb and Thirring |LT| . Lieb-Thirring inequalities refer to estimates 
that bound moments of negative eigenvalues of Schrodinger type operators in terms of the 
external fields. They play a fundamental role in various results concerning localized many- 
fermion systems. Most notably, the ground state energy of the many-body Hamiltonian 
in many cases is related to the sum of the negative eigenvalues of an effective one-body 
Hamiltonian. Among other useful applications, Lieb-Thirring inequalities stand behind the 
most effective and elegant proofs of stability of matter. They also serve as a powerful apriori 
estimate for the semiclassical analysis of the many-fermion ground state. 

We shall not attempt to give an overview of this vast and beautiful subject, since it is 
much better to refer the reader to the concise review article of Elliott Lieb [CTlj . We focus 
instead on the particular case of magnetic Lieb-Thirring (MLT) inequalities. They estimate 
moments of negative eigenvalues ei(H) < e.2(H) < . . . < of the Pauli operator 

H := [cr ■ (— zV + A)] 2 + V (1.1) 

on L 2 (R 3 , C 2 ) with a vector potential A, magnetic field B := V x A and external potential 
V. Here cr ■ v = a l v% + o~ 2 V2 + o~ 3 Vs, v G R 3 , and cr 1 ,^ 2 ,^ 3 are the Pauli matrices. Unlike 
in the nonmagnetic case, where the optimal form of the estimates is well-established and the 
remaining main challenge is to find the optimal constants, the magnetic case is much less 
understood. Apart from the case of the constant magnetic field, which has been settled in 
|LSY-llj (in two dimensions |LSY-lll| ). so far there is even no conjecture for a general magnetic 
Lieb-Thirring inequality that could be considered optimal in all aspects. 

Magnetic fields created in laboratories are usually weak and with a fair precision they can 
be handled perturbatively. Strong magnetic fields, however, occur in certain astrophysical 
models (e.g. neutron stars, see |LSY-1| ILSY-llj and references therein). Even a physically 
weak magnetic field can become effectively strong in certain problems related to quantum 
dots |LSY-lllj . 

Both these physical applications and the spirit of universality encompassed in the original 
Lieb-Thirring inequality have led us to search for magnetic Lieb-Thirring inequalities on the 
sum of the negative eigenvalues, that are optimal as far as the field strength is concerned in 
the strong field regime. 
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The semi classical formula for the sum of the negative eigenvalues, \ej(H)\, for a con- 
stant magnetic field behaves linearly in the field strength, |B| ( |LSY-IIj ). This fact suggests 
that Ylj \ e j{H) \ may be bounded by an expression that grows with the first power of |B| even 
for nonconstant magnetic fields and away from the semiclassical asymptotic regime. Our goal 
is to establish such MLT estimates with as few technical assumptions on B as possible and 
no technical assumptions on V. 

In this paper we present two such estimates. The simpler bound, Theorem 12.11 is proven 
in this paper. The more involved bound, Theorem 13.21 is outlined in Section but the details 
of the proof appear elsewhere jES-IV j. We point out that the methods behind these two 
proofs are very different and they are somewhat complementary, however, we did not succeed 
in combining the merits of both. 

While both bounds are optimal as far as the potential and the strength of the magnetic 
field are concerned, they require additional technical assumptions on the magnetic field. These 
are usually formulated in terms of variation lengthscales, and practically they are regularity 
assumptions on B. This means that supremum norms of derivatives of the magnetic field 
appear in the final Lieb-Thirring inequality. 

The difference between the two theorems is that the more involved bound, Theorem 13.21 
involves only local supremum norms. Therefore it enjoys an important locality property: the 
estimate is insensitive to the behavior of the magnetic field far away from the support of [V]_, 
where [a]_ := — min{0,a} denotes the negative part of a. The simpler bound, Theorem 12. 1[ 
involves the global C 5 -norm of the direction of the magnetic field, n := B/|B|. In particular, 
irregular behaviour of n far away from the support of [V]_ renders our estimate large despite 
that it should not substantially influence the negative spectrum. As a compensation, we need 
less assumptions on the regularity of the field strength |B| and the proof is much shorter. 

Finally we remark that armed with a MLT inequality that scales linearly in the field 
strength, one may prove the naturally expected semiclassical asymptotics for the sum of the 
negative eigenvalues of 

H(h, b) := [<x • (-ihV + bA)] 2 + V 

uniformly in b as h — > 0. This requires combining the techniques of ES-II and |ES-IVj and 
the details will be published separately. 

1.2 Short history of magnetic Lieb-Thirring inequalities 

For a constant magnetic field, B = const, the inequality 

\ej(H)\ < {canst) [J |B|[y] 3 _ /2 + jjvf! 2 ) (1.2) 
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proven in |LSY-llj is optimal, apart from the constant. It has seemed to be reasonable 
to conjecture that (J1.2)) also holds for an arbitrary magnetic field. However, such a naive 
generalization fails for two reasons. 

Firstly, even when B has constant direction in R 3 (jl.2|) can be correct only if |B(x)| 
is replaced by an effective field strength, B e s(x), obtained by averaging |B| locally on the 
magnetic lengthscale, |B|~ 1//2 ( E-1995J. 

Secondly, the existence of the celebrated Loss-Yau zero modes |LYj contradicts (jl.2j) . 
Indeed, for certain magnetic fields with nonconstant direction the Dirac operator T> : = 
(7 ■ (— zV + A) has a nontrivial L 2 -kernel. In this small potential perturbation of 

V 2 shows that ^ - \ej(H)\ behaves as f n(x)\V(x)]-dx, i.e. it is linear in [V]-. Here n(x) is 
the density of zero modes, n(x) = 52 . \ u ji x )\ 2 i where {uj} is an orthonormal basis in Ker£>. 
Thus an extra term linear in [V]_ must be added to (jl.2|) and n(x) has to be estimated. 

The problem of the effective field, observed in [E-1995J, was first succesfully addressed by 
Sobolev, [5ob-96j. [S ob-97j and later by Bugliaro et. al. [BFFGH] and Shen |Shj . In particular, 
the L 2 -norm of the chosen effective field, ||-B e ff|| 2 , is comparable to ||B|| 2 in BFFGS, and the 
same holds for any L p -norm in Shen's work. In a very general bound proved in |LLSj the first 
term in (jl.2|) is replaced with ||B||2^ 2 || V^|| 4 . 

In the works jE-1995| . |Sob-97| . |BEj . [EES] . |BFFGS| the density n(x) is estimated by a 
function that behaves quantitatively as |B(x)| 3 / 2 . In particular, in the strong field regime 
these estimates are not sufficient to prove semiclassical asymptotics for H(h, b) uniformly in 
b; they give the asymptotics only up to b < (const. )h -1 Sob-98j. 

We remark that the bounds in |LLSj and [BFFGS have nevertheless been very useful in 
the proof of magnetic stability of matter. In this case the magnetic energy, J |B| 2 , is also part 
of the total energy to be minimized, therefore even the second moment of the magnetic field 
is controlled. We also remark that if the field has a constant direction, then no Loss-Yau zero 
modes exist, n(x) = 0. In this case Lieb-Thirring type bounds that grow linearly with |B| 
have been proved in |E-1995j and |Sob-96j . |Sob-97j . 

Therefore the fundamental difficulty is to understand the density of Loss-Yau zero modes. 
It is amusing to note that it was a substantial endeavour to show that a zero mode may exist 
at all |LYj , and that multiple zero modes may also occur jES-IIIj . On the other hand, it seems 
also quite difficult to give an upper bound on their number in terms of the first power of the 
field strength ( Corollaries 12 . 21 and 13 . 3j) . 

Since n(x) scales like (length) -3 and |B(x)| scales like (length) -2 , a simple dimension 
counting shows that n(x), therefore Ylj l e j|> cannot be estimated in general by |B(x)| or by 
its smoothed version. However, if an extra lengthscale is introduced, for example certain 
derivatives of the field are allowed in the estimate, then it is possible to give a bound on the 
eigenvalue sum that grows slower than | !B | 3 / 2 in the large field regime. Before the current 
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paper and our recent work |ES-1V| there were only two results in this direction. 

The work |BFGj used a lengthscale on which B changes. The estimate eventually behaves 
like |B| 17 / 12 in the field strength. As far as local regularity is concerned, only B 6 H\ oc is 
required. However, n{x) is estimated by a quantity that depends globally on B(x) not just in 
a neighborhood of x. As a consequence, no locality property holds for the MLT inequality in 

[EES]. 

Our earlier work |ES-lj had a different approach to reduce the power 3/2 of |B| in the 
estimate of n(x). We introduced two global lengthscales, L and £ respectively, to measure the 
variation scale of the field strength |B| and the unit vector n := B/|B| that determines the 
geometry of the field lines. This required somewhat more regularity on B than [BFG*j and 
it also involved the W l,1 -norm of V. The estimate grew with the 5/4-th power of the field 
strength |B| in the large field regime. For fields with a nearly constant direction, i 3> 1, the 
bound was actually better, it behaved like |B| + |B| 5//4 £ -1 / 2 . This indicates that it is only 
the variation of n and not that of |B| that is responsible for the higher |B|-power. The MLT 
estimate in jES-Ij did not enjoy the locality property either. 

Due to the improvement in the |B|-power from 3/2 to 5/4 in the MLT estimate we could 
also prove the semiclassical eigenvalue asymptotics for H(h, b) in the regime b <C h~ 3 for 
potentials in W 1 ' 1 |ES-11| . This bound turned out to be sufficient to show that the Magnetic 
Thomas- Fermi theory exactly reproduces the ground state energy of a large atom with nuclear 
charge Z in the semiclassical regime, i.e. where b -C Z 3 , Z oo ES-II]. The condition 6<C 
is optimal as far as the semiclassical theory is applicable. 

Despite the successful application of the bound in [BFGj to the stability of matter with 
quantized electromagnetic field with an ultraviolet cutoff [BFrTj] . and despite that the MLT 
inequality given in jES-Ij fully covered the semiclassical regime of large atoms, it is still 
important to establish a uniform Lieb-Thirring type bound with the correct power in the 
magnetic field and with no unnecessary assumptions on V. Such bound will likely be the 
key to generalize the analysis of the super-strong field regime of jLSY-Il to non-homogeneous 
magnetic fields. 



2 Lieb-Thirring inequality without locality property 

We consider the three dimensional Pauli operator, H = V 2 + V, with a differentiable magnetic 
field B = V x A. The operator 

V := a ■ (— iV + A) (2.1) 

is called the Dirac operator with vector potential A and magnetic field B. We make two 
global assumptions: 
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Assumption 1. B(x) 7^ for all x G R 3 , i.e. the unit vectorfield n := B/|B| is well defined. 
Assumption 2. The vectorfield n satisfies the following global regularity condition 

5 

L' 1 := E HV 7 n||^ < 00 , (2.2) 

7=1 

where L n is called the global variation lengthscale of n. 
For any L > 0, x G R 3 we also define 

B* L (x) :=sup{|B(y)| : |y - x| < L} + L ■ sup{| VB(y)| : |y - z| < L} . 

Theorem 2.1 (Lieb-Thirring inequality without a locality property). For any < 

L < -^n; swm 0/ the negative eigenvalues, ei(H) < 62(H) < . . . < 0, o/iif satisfies 

J2 \ej{H)\ < c^L- 1 J (Bl + L- 2 )[V]- + 1 52[V] 3 _ /2 + J [V] 5 / 2 ) (2.3) 
a universal constant c. 

Let Ht> 2 <\ denote the spectral projection onto energy levels below A in the spectrum of 
T> 2 . Along the proof of Theorem 12.11 we obtain the following bound on the density of the low 
lying states of V 2 : 

Corollary 2.2 (Density of low lying states). For any < L < L n and x G R 3 

np2< A (i, 1) < c\ l J 2 m&x{B1(x), A*} , with A* := max{A, L~ 2 }, (2.4) 

in particular the local density of zero modes of T> 2 grows with at most the first power of the 
field strength. 

Convention. We use the letter c to denote various positive universal constants whose exact 
value may change from line to line. 

3 Lieb-Thirring inequality with a locality property 

In this Section we give an outline of a more refined magnetic Lieb-Thirring inequality. The 
detailed proof appears elsewhere lES-IVj . 



We assume that B G C 4 (R 3 ,R 3 ), and define three basic lengthscales of B. The Pauli 
operator will be localized on these lengthscales. Let n := B/|B| be the unit vectorfield in 
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the direction of the magnetic field at all points where B does not vanish. For any L > and 
x G R 3 we define 

B L (x) :=sup{|B(y)| : \x - y\ < L} , (3.1) 

and 

b L (x) :=inf{|B(j/)| : \x - y\ < L} (3.2) 

to be the supremum and the infimum of the magnetic field strength on the ball of radius L 
about x. 

Definition 3.1 (Lengthscales of a magnetic field). The magnetic lengthscale o/B is 

defined as 

L m (x) := sup{L > : B L (x) < L~ 2 } . 
The variation lengthscale o/B at x is given by 

L v (x) := min{L s (x),L n (x)} , 

where 

L s (x) :=su P {l>0 : LP sup { | V 7 |B(y)| | : \x - y\ < i,} < b L (x), 7 = 1, 2, 3, 4} 

L n (x) :=sup{L>0 : I? sup{|Wn(j/)| : |x - y| < L, B^) ^ 0} < 1, 7 = 1, 2, 3, 4} 
fwf/i i/ie convention that sup0 = — 00). Finally we set 

L c (x) := max{L m (x) , L v (x)} . (3.3) 

A magnetic field B : R 3 — > R 3 determines three local lengthscales. The magnetic length- 
scale, L m , is comparable with |B| -1 / 2 . The lengthscale L s determines the scale on which the 
strength of the field varies, i.e., it is the variation scale of log |B|. The field line structure, 
determined by n, varies on the scale of L n . The variation lengthscale L v is the smaller of these 
last two scales, i.e., it is the scale of variation of the vectorfield B. 

For weak magnetic fields the magnetic effects can be neglected in our final eigenvalue esti- 
mate, so the variational lengthscale becomes irrelevant. This idea is reflected in the definition 
of L c ; we will not need to localize on scales shorter than the magnetic scale L m . 

The following theorem is the main result in [ES-IVJ. 

Theorem 3.2 (Lieb-Thirring inequality with a locality property). We assume that 
the magnetic field B = V x A is in C 4 (R 3 , R 3 ). Then the sum of the negative eigenvalues of 
H = [er ■ (-iV + A)] 2 + V satisfies 

5>;|<c|[F] 5 _/ 2 + c J \B\lvf 2 + cJ(\B\ + L-i)L- l [V}_. (3.4) 
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Corollary 3.3 (Density of zero modes with a locality property). Given a magnetic 
field B e C 4 (R 3 ,R 3 ) ; the density of zero modes of the free Dirac operator T> with magnetic 
field B satisfies 

n(x) := \uj(x)\ 2 < c(\B(x)\ + L;\x))L;\x) , (3.5) 
j 

where {uj} is an orthonormal basis in the kernel ofV = cr- (— zV + A). 

Remarks, (i) The bound ()3.5|) is optimal as far as the strength of the field |B| is concerned. 
This fact follows from the construction of Dirac operators with kernels of high multiplicity 
following the method of F.S-fTT]. 

(ii) Notice that the Lieb-Thirring inequality of |LSY-IIj for a constant field is recovered 



in Theorem 13.21 However, the uniform Lieb-Thirring bound for a constant direction field, 
|Sob-97j . |ES-Ij . does not directly follow from our main theorem as it is stated. Firstly, ()3.4|) 
contains a term linear in V that is unnecessary for a constant direction field. Secondly, we 
assume high regularity on B. This regularity is needed only to construct the appropriate 
curvilinear cylindrical localization, which is unnecessary for a field with constant direction. 

(iii) It can be shown that L c (x) enjoys a locality property; for any 5 > the inverse 
lengthscale L~ 1 (x) is bounded by a ^-dependent function of the local C 4 -norm of n and 
log |B| in a 5- neighborhood of x. 

Before we turn to the proof of our new result, Theorem 12.11 we explain the key difference 
between the two proofs that is somewhat hidden behind technicalities. 

The linear power of |B| in the estimate reflects the basic fact that the space with a magnetic 
field cannot be considered isotropic: the quantum motion parallel with the magnetic field 
behaves differently than the transversal one. The magnetic field affects only the transversal 
motion; this is why the eigenvalue sum scales in the same (linear) power of |B| both in d — 2 
and 3 dimensions. The motion in the third direction affects only the powers of the potential, 
as it is customary in the nonmagnetic Lieb-Thirring inequalities. 

All MLT estimates that yield |B| 3 / 2 behaviour neglect this geometric fact by simply com- 
paring the magnetic problem with a nonmagnetic one, usually via a diamagnetic inequality 
that loses the anisotropic feature of the problem. The typical estimate is of the form 

V 2 > b -i V 2 = 6-i(D 2 + o--B) (3.6) 

where D := — z'V + A is the spinless magnetic momentum and b := ||B|| 3> 1 denote some 
(local) norm of B. The Pauli kinetic energy is scaled down so that the dangereous cr ■ B 
becomes bounded uniformly in b. The magnetic Laplacian D 2 then can be controlled by the 
nonmagnetic Laplacian, —A, but the factor b^ 1 now affects all three coordinates, yielding a 
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scaling of b 3 ^ 2 . The key to the improvement is to separate the motion in parallel and in the 
transversal directions and use a crude estimate similar to (|3.6|) only in the two-dimensional 
transversal kinetic energy. 

Since the direction of the magnetic field n varies, it is in general impossible to introduce 
global orthogonal coordinates parallel and transversal to n. The more conventional approach 
is the one presented in [ES-IV : we introduce appropriate coordinates locally and approximate 
the magnetic field by another one that is constant in these coordinates. For constant magnetic 
field the separation is then trivial. The main technical challenge is to keep all approximation 
and localization errors bounded uniformly in the field strength b. This requires many steps: 
most importantly we have to change the metric by a conformal factor to make the field strength 
constant along the field lines; we have to localize the problem onto the physically correct 
domains, which are elongated curvilinear cylinders of width 5 -1 / 2 along the field lines and 
we have to construct appropriate coordinate systems along each field line. The construction 
steps require tools from the spinor geometry behind the Dirac operators. The localization 
onto narrow cylinders is done via the magnetic localization formula from jES-IIj . All these 
steps require fairly high regularity on the magnetic field. 

The more direct approach is the one behind Theorem 12.11 This method bypasses all the 
complications with the constant field approximation, the cylindrical construction and the 
change of metric. It uses the simple observation that the energy in the parallel direction can 
be extracted from T> 2 by a simple energy argument (see (j6.18|) in Lemma f6.2|) : 

V 2 > cD 2 n - (error) , (3.7) 

where _D n := — id n — |divn, and the error terms are of lower order in momentum. It turns out 
to be necessary to consider higher powers of the resolvent (T> 2 + A) -1 to control the ultraviolet 
regime (here A > is a constant controlled by the geometry of the field lines), therefore it 
is necessary to estimate higher powers of T> 2 as well, the key observation is that after having 
extracted D^, higher powers of the crude diamagnetic estimate (|3.6|) are sufficient. Although 
()3.6|) cannot be squared directly, after several commutators it is possible to show that (see 

(JSD31) 

where S := min{A6 _1 , 1}. 

Therefore the basic estimate for the density of the low lying states is 

8A 3 cA 3 cA 3 

lll)2<A < , \ \ q < , ao^o , \ q < 



(V 2 + A) 3 ~ V e + \ 2 Dl + A 3 - (5D 2 ) 3 + XWl + A 3 ' 

Then we can essentially use the diamagnetic inequality (although first we have to ensure that 
the two operators in the denominator commute, see Lemma 16.1)1 and the diagonal term is 
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estimated as 

It is then easy to change this quantity into Euclidean coordinates and compute 

A 3 



/ X3f 2^3 , X2 2 , X3 d P ^ c\ 1/2 (b+X) 

7r3 5 3 (p 2 ) 3 + A 2 ^ + A 3 



which scales linearly with b. 

Apart from several technicalities, there is an apparently innocent but fundamental compli- 
cation in both approaches. Both proofs uses the concept of the local magnetic field strength, 
b, that is given by some local norm of B. Therefore at the beginning one needs to localize the 
original operator T> 2 onto domains of size independent of b to establish this apriori localization 
scale. We choose domains of size L n where the field lines do not vary substantially. 

Localizing higher powers of T> 2 even onto the scale L n is difficult because it does not follow 
simply from energy considerations. If the localization length L n is a uniform constant then it 
can be done using the spectral theorem (see (|4.2j) ). This is the method in Theorem 12.11 but 
it forces one to introduce the same localization error everywhere in space, even in regimes, 
where the field is very regular. 

If one wants to include a variable L n , as in Theorem 13.21 to ensure the locality property, 
then one needs a more powerful apriori localization. Typically it is not hard to localize 
resolvents of second order elliptic operators onto cubes of size i at the expense of an error 
£~ 2 . However localizing the square (or higher powers) of the resolvent requires off-diagonal 
estimates on the resolvent kernel (see Proposition 7.1 in ES-IV ). While these are typically 



easily available for scalar elliptic operators without spin, we do not know any apriori off- 
diagonal control on the resolvent of V 2 . If the original Pauli operator is estimated by a 
constant field Pauli operator, then aposteriori we can extract off-diagonal estimates, and this 
idea is used, though very implicitly, in |ES-IVj . But without comparison with the constant 
field problem, we do not have off-diagonal control. This is the main reason why we are unable 
to extend the elegant and short method of Theorem 12.11 to give any locality properties. 



4 Organization of the proof of Theorem 12.11 

We can assume that the potential in ()1.1|) is nonpositive, so for simplicity we can consider 
H = T> 2 — V with V > 0. We start with the Birman-Schwinger principle 

n (\ v - e \- 2 ^Te iv ~ El ~ 2 A ) dE > (41) 
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where n(X; c) denotes the number of eigenvalues of the positive self-adjoint operator X above 
level c. 

We estimate T> 2 from below by the magnetic Laplacian, then use the diamagnetic inequality 
to estimate its resolvent kernel by that of the free Laplacian. These estimates involve error 
terms that cannot be controlled by E in the resolvent (T> 2 + E)" 1 when E is small. As a first 
step, we insert a positive constant P in the resolvent using the inequality 

1 2 8 P 3 . , 

< ^T-TT-^ + -B 7^> ■ m, (4-2) 



V 2 + E ~ V 2 + P + E E(V 2 + Pf 

that follows from the spectral theorem and the corresponding arithmetic inequality. P is 
chosen as 

P := e~ 4 L- 2 , 

where e denotes a sufficiently small universal constant. 

From ()4.2j) and the general bound n(Xi + X 2 ; c\ + c 2 ) < n{X^ c\) + n(X 2 , c 2 ), we obtain 



with 



\Tr(V 2 - V)_| < c[(J) + (II)] (4.3) 



/ /I \ 

^-E^^^W-Etf.^iE; (4.4) 



(II) -=1 n(\V- El^jp^ylV ~ E] 1 / 2 , ±E)dE . (4.5) 
These two terms will be called positive energy regime and zero mode regime, respectively (see 



ES-IV]). 



The proof of Theorem 12.11 then follows from the two estimates below: 
Proposition 4.1. For sufficiently small e 

(/) <c J Bl(x)V 3/2 (x)dx + c J V 5/2 (4.6) 

and 

(II) < cP 1/2 J msx{Bl(x), P}V(x)dx . (4.7) 

To obtain these estimates, we will have to bound the resolvent of T> 2 and the diagonal 
kernel of the cube of the resolvent of T> 2 (see (|7.2|) ). In both cases we localize the operator T> 2 
onto cubes of size of order L^ 1 . 
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In Section |S] we define localization functions and new Dirac operators which coincide with 
T> locally but have a uniformly bounded magnetic field everywhere. In Section El we prove 
various estimates on a Dirac operator with a bounded magnetic field and we compare it with 
an auxiliary Laplacian-type elliptic operator that is diagonal in the spin space. Finally, in 
Sections [7| and |H1 we complete the estimates ()4.7|) and (|4.6j) . respectively. 

5 Localization 

Consider the cubic lattice A = eLZ 3 . For each fee A, let Q k C fi' fe C Q* k C fif be the balls 
with center k and with radius eL, 2eL, 3eL and 4eL, respectively. Let {6k}keA be a partition 
of unity such that < 9 k < 1, 

J2 d l = ^ supple fi fe , ||V^fe|U < c{eL)-\ 1<7<5. (5.1) 

fceA 

We also choose further localization functions i] k , uj k and Xk with 9k < % < ^fe < Xfc < 1, 
A; G A, such that 

supp 7] k cQ' k , r] k 

supp u; fc C fife, ^ 
supp Xfc C fif , Xfc 
We note that 

£><C, $> fc <C, J>*<C (5.5) 

fceA fceA fcGA 

for some universal constant C by the finite overlap property of the balls. 
Let b k := sup n #(|B| + L|VB|), then we have 

k 

|B| + L|VB| < b k < B* L on supp(x fc ) , (5.6) 

if e is sufficiently small. 

The following lemma defines a new Dirac operator for each cube Q k that coincides with D 
on fi£ but has a uniformly bounded magnetic field. The proof is given in Section IA.1I 

Lemma 5.1. For each k G A there exists a nowhere vanishing magnetic field B k such that 
B fc = B on Q* k , Bfe = (const) on R 3 \ flf, 

HBfelU + LUVBfelU <cfefe, (5.7) 



= 1 


on 


fifc, 


l|V 7 r/ fc || 


= 1 


on 




||V^ fc | 


= 1 


on 




W^Xk 



1 < 7 < 5 ; 
1 < 7 < 5 ; 



(5.2) 

(5.3) 
(5.4) 
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and rifc := B^/|Bfc| satisfies the bound 

5 

£(eL)lV 7 n fc || 00 <ce. (5.8) 

7=1 

Moreover, there exists a vector potential A^, V x = such that = A on Vt* k . The 
Dirac operator T>k '■= cr ■ (— iV + A&) magnetic field therefore satisfies 

l^kVk = Vri k . (5.9) 

The localization is performed by the following "Pull-up" formula whose proof was given 
in BFFGS and was also recalled in Proposition 6.1. in |ES-IVj . 

Proposition 5.2 (Pull-up formula). Let I be a countable index set and let gi, i G / , be a 

family of nonnegative functions on R 3 such that Yliei 9i( x ) < 00 f or every x G R 3 . Let Ai, 
i G /, be a family of positive invertible self-adjoint operators on L 2 (R 3 ; C 2 ). Then 

(E^E^W(^ 9 £ E^ft- I"") 
The resolvent (V 2 + P + E)- 1 can be easily localized using the partition of unity {6k}- 

fc fc k k 



using (|5.9)l . the finite overlap property of the supports of 9k s, and the estimate on their 
derivatives (|5.1j) . For sufficiently small e the error can be absorbed into P, so 



V 2 + P + E > ]^9 k {V 2 k + P + E)9, 

k 

Using the "Pull-up" formula (J5.1U)) . we obtain 
1 2 

< = : < 



V 2 + P + E -= J2k8k(Vl + P + E)6 k " 2 ^ 6k Vl + P + E 9h - (511) 

Similar localization does not hold for higher powers of the resolvent as it was remarked in 
Section 7 of [ES-IV . However, we can localize the second term in ()4.2|) estimating it by an 
auxiliarty elliptic operator that is independent of the spin coordinates. The construction is 
given in Section El 
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6 Analysis of bounded magnetic fields 



Throughout this section we assume that a vector potential A is given so that for some positive 
constants L, e and b, the nowhere vanishing magnetic field B = V x A of the Dirac operator 
T> = a ■ (— zV + A) satisfies the following conditions: 

B = const outside of (6.1) 

where f2* is a ball of size 0(eL); 

< |B| +L|VB| < b (6.2) 

and 

5 

J2(^y\\Vn\\oo < ce , (6.3) 

7=1 

where n := B/|B|. The plane V through the center of f2# and orthogonal to the magnetic 
field at the center is called the supporting plane. 

The results of this section will be applied to the Dirac operators T>k with a magnetic field 
Bfc constructed in Lemma f5. 11 but we do not carry the index k inside this section. 

We choose a global, positively oriented orthonormal basis {e\,e2,e%\ in R 3 whose third 
vector is e 3 := n, and such that all Christoffel symbols T 7 ^ = (V e ^e&, e m ) vanish outside the 
ball twice bigger than and they satisfy the estimate 

4 

^Ln^^Woo < cL- 1 . (6.4) 

7=0 

That such and othonormal basis exists follows easily from (|6.3jh e.g., by a Gram-Schmidt 
procedure starting from n, ei,e2, where ei,e 2 are fixed vectors. We will work on the trivial 
spinorbundle over R 3 , i.e., in the L 2 (R 3 ) ® C 2 space. We briefly recall the appropriate 
formalism from Section 9 of [ES-IV . 

For any 1-form a we define the spinor connection 

V£ := d x + ia{X) + itr ■ uj{X) , 

where 

ui{X) := ^(Vxe 3 ,e 2 ), (Vxei,e 3 ), (Vxe 2 ,ei)j , 

and let 

II?:=-iV£, J = 1,2,3. 
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For simplicity we will use Vj := V e and dj = d e . and we usually omit the superscript a from 
the notation. We recall that the operator 

D x := -id x - f div X (6.5) 

is self-adjoint on L 2 (R 3 ) for any C 1 bounded vectorfield X, therefore the operator 

D x := -id x - f div X + a(X) (6.6) 

is self-adjoint. The hat refers to the operators without magnetic field, a = 0. 

These operators are originally defined on L 2 (R 3 ), but with a little abuse of notations they 
can also be viewed as acting on L 2 (R 3 ) ® C 2 as Dx ® I2, where I2 is the identity operator on 
the spin space. We will not distinguish between these two operators in the notations. We set 
Dj := D e . for simplicity, then 

IT, = Dj + \ div ej + \a ■ u( ej ) . (6.7) 

We compute the commutators of the -D/s: 

\D^D k } = {-t){D M +(3 3k ) (6.8) 

where (3 := da and [3jk '■— /3(ej,ek). In the computation we used that the divergence part 
vanishes since 

Xdiv Y - Ydiv X - div [X, Y\ = (6.9) 

for any vectorfields. The relation (J6.9)) is easy to check on coordinate vectorfields and it is 
tensorial. 

Let Qjl := r™j — Tjl, then D[ e . >eh ] = \{Q^ k D m + D m Q^.), where we used the summation 
convention for repeated indices. From (|6.8J) 

[Dj, D k ] = -l(Qf k D m + D m Qf k ) - i(3 jk = (-i) (q? k D m - l -(d m Qf h ) + (3 jk ) . (6.10) 

We recall from Section 9 of |ES-IVj that is a spinor connection with a magnetic 2-form 
(3 := da. If we choose the vector potential a(X) := A-X, then ~kf3 is the 1-form corresponding 
to the vectorfield B, where * denotes the Hodge dual with respect to the standard Euclidean 
metric. With these notations the Dirac operator T> = cr ■ (— zV + A) in Euclidean coordinates 
is SU (2)-gauge equivalent to 

V:=a-n, (6.11) 

i.e., 

V = UDU* (6.12) 
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for some U : R 3 i— > SU(2). Tilde always refers to Dirac operators written in a basis where 
the third basis element is parallel with the magnetic field. In the chosen basis the magnetic 
2-form satisfies /3(ei, e^) = B, where B :— |B| is the field strength, and all other /3(ej, ej.) = 0. 
The Lichnerowicz formula is given by 

v 2 = n* • n + -Rq + o-(*/?) , (6.13) 

where Ro is the scalar curvature (which vanishes in our case) and for any one-form A = Ylj 
we define <r(A) = J2j a ^j- 

Let Z be a symmetric positive definite 3 by 3 matrix valued function on R 3 that satisfies 

3 

\\Z-I\\ +^(£L) 7 ||V 7 Z|| <ce, (6.14) 

7=1 

and let h be a function with 

2 

^(^llV^Hoc < ce^L' 2 . (6.15) 

7=0 

We define 

T := T ZA := D -ZB + h, f := f ZA := D • ZD + h (6.16) 

on L 2 (R 3 ), but we will also denote simply by T, T the operators T ® h and T ® I 2 acting on 
L 2 (R 3 ) ® C 2 . We will choose Z and h according to the following lemma whose proof is given 
in Section fA. 21 

Lemma 6.1. Let the magnetic field B = V x A satisfy \6.1\) . \6. ty) . and \6. ty) . Then for 
sufficiently small e there exists a 3 by 3 symmetric matrix valued function Z and a real 
function h satisfying \6. 14\ ) and h6.1fy) . such that the commutator of T = D • ZT) + h and D n 

vanishes, 

[T,D a ) = 0. (6.17) 

Armed with these notations, we state the following lower bounds on the Pauli operator. 
The proof is given in Section IA.31 

Lemma 6.2. Let the magnetic field B of the Dirac operator D = cr ll satisfy \6.1\) . \6. jjj) . and 

JO) . Let < 9 < 7] < u < x < 1 be real functions on R 3 with |V 7 6»|, |V 7 r]|, |V 7 u;|, |V 7 x| < 
c(eL)~ 7 ; 1 < 7 < 5 ; and rj = 1 on supp(6), uj = 1 on supp(n) and x = 1 on supp{uj). Let 
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P := e 4 L 2 , A > P and 5 := min{A6 l , 1}. Then for sufficiently small e there exists a 
positive universal constant c such that 

V 2 >cD 2 n -e 2 P (6.18) 

and 

V 3 6 2 V 3 + XV 2 r) 2 V 2 + X 2 Voo 2 V + X 3 X 2 > c6(5T) 3 6 , (6.19) 
where T = Tz,h is defined in \6.1b]) with the bounds \(j-14\l , \b'-lty - 



Finally, we estimate the diagonal of the resolvent kernel of the sum of the lower bounds 
obtained in Lemma 16.21 The proof is given in IA.4I 



Lemma 6.3. Let the magnetic field B = V x A satisfy \6.1\) . \6. ty) . and \6. and let 

A > e~ 4 L~ 2 , 5 := min{A6 _1 , 1}. Then for sufficiently small e 

x, x) < cA 1/2 max{6, A} . (6.20) 



(5T) 3 + X 2 D 2 n + A 3 

7 Proof of the zero mode regime 

In this section we prove a Lemma that immediately implies (J4.7|) . The main tool is a local- 
ization of the third power of the resolvent. For later purposes, we will need the following, 
somewhat more general estimate: 

Lemma 7.1. Assume that the magnetic field B of a Dirac operator T> = a ■ (— zV + A) is 

nowhere vanishing and that it has a finite, nonzero global variation lengthscale L n Let 
L < L n , then for sufficiently small e and for any X > P = e~ A L~ 2 we have 

A 3 

'x,x) < cX 1/2 m&x{B* L (x),X}. (7.1) 



(V 2 + A) 



For the proof of ()4.7|) we can use n(XY 2 X; c) = n(YX 2 Y; c) for positive operators X, Y, 
to bring \V — E\_ in the middle in (J4.5|) . and estimate it by V, then undo it 

UD < f 4 vll2 (v^py v1 ' 2 • ^ E ¥ E = 16P ' Tl ( y1 ' 2 why y1 ' 2 ) ■ (7 - 2) 

Using (|7.1|) for A := P, we immediately obtain ()4.7|) . Using that U^ 2<x — SX 3 (V 2 + A)~ 3 from 
the spectral theorem, we also obtain Corollary 12.21 Note that ()2.4|) holds for any A > since 
the case A < A* is trivial. 
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Proof of Lemma \7.1\ We insert the partition of unity, recall the bounds (|5.5|) and use (|5.9|) 
to obtain 



(V 2 + A) 3 = [^ 3 ^^ 3 + 3\V 2 9 2 k V 2 + 3X 2 V9 2 k V + A 3 ^ 

fceA 

> c [VktiPk + M%n»Pk + ^kC0 2 k V k + \ 2 V k 9 2 k V k + A 3 ^] • (7.3) 
fceA 

i 2 TV /)2 



We use the X T> k 9 k T> k term to save a derivative in the third direction. Let T> k be the Dirac 
operator constructed in SectionElfor the vector potential A k . Then there exists a spin rotation 
U k : R 3 i-f 5C/(2) (see lj?H2|» ) such that 

P fe = U k V k U; . (7.4) 

By Schwarz' inequality, |Vfc#| < c(£L) _1 and (j6.18j) from Lemma f6.2l applied to the oper- 
ator T>\ we have 

V k 9 2 V k > c9 k V 2 9 k - c(eL)- 2 Vk > c9 k U k D 2 nk U* k 9 k - ce 2 P Vk . 

We can sum up these inequalities, use (|5.5|) and A > P, to obtain from (|7.3|) that for sufficiently 
small e 

(t? 2 + a) 3 > cJ2[vleM + \v 2 kV 2 k v 2 k + \ 2 v k u 2 k v k + \ 3 xl 

k 

+cX 2 J2 e kU k (D 2 nk + X)U* k 9 k . 
k 

Applying the estimate (j6.19J) from Proposition 16.21 to the operators V k = U k V k U k we obtain 
that 

(X? 2 + A) 3 > d kU h ((4T fe ) 3 + \ 2 Dl k + A 3 ) U* k 9 k , 

k 

where 5 k := minjAfo^ 1 , 1}, the operators D k and T k are obtained by applying the construction 
of Section IH1 for the vector potential A k . 

Hence by the "Pull-up" formula ()5.1U|) we get 

1 S c £ e k U„ ^ >3 DJfc . (7.5) 



We apply ()6.20|) for the magnetic field B fc = V x A fc satisfying the bound ()6.2|) with b = cb k 
from ()5.7|) to get 

— — — < cA 1/2 max{6 fc ,A}, 
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and use that only finitely many 9 k s are nonzero at any x e R 3 , and for all such fc's we have 
b k < B* L (x). This completes the proof of (J7.1j) . □ 

8 Proof of the positive energy regime 

In this section we prove (|4.fi)) . We start by applying the localization (|5.11jl in (|4.4jl 

(/) < [ »(e i" - ^^rrk^ 11 ' - E| - 2 • ! 

Since the supports of at most C of the functions 6 k can overlap (see (|5.5Jl ). we can "pull out" 
the summation at the expense of descreasing the energy threshold \ (see Section 6 of |ES-IVj 
for more details). Therefore 

n(\V - E\]!\ w -—-— 9 k \V - El 1 / 2 , j^)dE . (8.1) 

Kfc/V - k 

For any A € R let 

II*(A) := X (^<A) 

denote the spectral projection below the level A in the spectrum of T>\. We have 

(I) < (J/1) + (J/2) 

with 

(//!):=£ r „(|V - g|l V + P " - , ^)<1E , (8.2) 

(//2):=£ /"n(|V - Etrth ^g-f-^ W - E\^ , ^)iE . (8.3) 

The estimate of the first term is easy: 

/ - U k (b k -P-E) < c < c 

V 2 k + P + E ~ V 2 k + cb k + E ~ (-iV + A fc ) 2 + £ 

using the spectral theorem first, then the Lichnerowicz formula in the Euclidean coordinates: 

v l = HV + A k ) 2 + tr ■ B k 
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and cr • Bfc > —cbk, where the constant is from (|5.7|) . Then we can proceed as in the proof of 
the original Lieb-Thirring inequality jLTj for each summand in (J8.2)) : 



n[\V - E\ l ! 2 e k 



(-iV + A k ) 2 + E 
< 64(C + 1) 2 



e k \v-E 



1/2 



1 



8(C+1) 



Tr( \V -E\ l l 2 e k 



dE 

c 





< c I Tr[ \V-E\ 2 _9 4 



-iV + A fc ) 2 + £ 

I 12 



1/2 



dE 



< c 



< c / \/ 5 / 2 



2 fl 4, 



|y(x)-£;| 2 _^ 



-iV + A k ) 2 + E 

1 12 



dE 



-A + El 



x, x)dxdE 



where we also used the diamagnetic inequality for the square of the resolvent. After the 
summation over k and using ([5.5)1 . we obtain the second term in the estimate ([4.6)1 . 

For the term (1/2), we first notice that the dE integration can be restricted to E < b k , 
and we can also assume that that P < b k - After this restriction we estimate 11^(6^ — P — E) 
by Hk(b k ) and we obtain 

- u k (b k ) a IT/ p|1/2 



(1/2) < [ k n(\V -E\ in 6 
fceA ^° 



k V 2 + P + E 



< S(C + 1)J2 Tr(\V - E^O, 

r.^A JO ^k 

" X(^<bk 



e k \v-E 
TikiPk) 



i 

' 8(C+1) 



dE 



fceA ' 



k V 2 + P + E 



9kW-E 



1/2 



dE 



fceA 



JR 3 



|V(a;) - E\_6l{x) 



X + P + E 



dU k (X) (x,x)dxdE , (8.4) 



using the spectral resolution of V 2 . Clearly 

X(A < h) 



X + P + E 



dU k (X) 



1 



< 



bk + P + E 
b^E UM + 



n fc (fe fc ) - 



i 



11,(0) + 



P + E 

P MP) 



(X + P + E) 2 

1 n fc (A)dA . 



n fc (A)dA 



p (A + E) 



•5) 



(P + E) 2 

The conclusion of Corollary 12 . 21 also applies to replacing the variation lengthscale L n with 
eL n , see ([5.8)1 . therefore 

U k (X)(x,x) < cA 1/2 max{S2(x), A} (8.6) 



20 



for A > e~ 6 L~ 2 , but then the same inequality also holds for all A > P = e~ 4 L~ 2 if c is changed 
to ce~ 3 using the monotonicity of A ^ n fc (A). 

From ([8.5)1 . ()8.6|) and the fact that P < bk < B* L (x) for x G supp (9k), we obtain 

< cs-*E- l l 2 Bl(x) . 
Therefore we can complete the estimate of (1/2) from (|8.4j) 

(1/2) < ce~ 3 V /°° / - E|_E- 1/2 dx d£ = ce~ 3 [ B* L V 3/2 . 

This completes the proof of ()4.6|) after choosing e to be a sufficiently small universal constant. □ 



A Proof of the technical lemmas 
A.l Proof of Lemma 15.11 

Let Vk be the plane through k £ A and orthogonal to n(k). We define := Xk n +(^~ Xk) n (k)- 
Since — n(k) = (n — n(k))xk — 0(e) from ()2.2|) . we have • n(A;) > 1 — 0(e) and |n^| = 
l + 0(e), in particular ^ for sufficiently small e. Moreover Vn^. = x^Vn+ (n — n(/c))Vxfcj 
therefore Vn^ = O(L^) since |n — n(k)\ = 0(e) on the support of \k from ()2.2|1 and L < L n . 
Similar bounds hold for the higher derivatives. We define := n^/ln^l and (J5.8|) can be 
easily checked. 

Now B fc is defined as B fc := f k n k where the function f k is constructed as follows. We 
define 

/*:=X*|B| + (1-X*)|B(A;)| Q nV k , 

then clearly \fk\ < cbk on Vk- 

Since the vectorfield is nearly parallel (clearly • n(A;) > 1 — 0(e)), its integral curves 
define a transversal foliation to the plane Vk- Hence we can extend fk from Vk onto the whole 
R 3 by ensuring that is divergence free, i.e., div B^ = fkdiv + V nfe /fc = 0. This requires 
integrating the equation 

V nfc (log/ fc ) = -div n fc 

along the integral curves of n fc . Since |Vn fc | < cL~ x , and it vanishes outside of Qf, we obtain 
\fk\ + L\V(fkiik) < cbk everywhere as required in ([5.70 . 
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The original field strength |B| also satisfies the equation 

V„(log|B|) = -div n, 

and since = n on Qk by definition and |B| = on 7^ H we obtain that fk = |B| on 
Q k . Therefore B fc = B on Q* k . 

This completes the definition of the extension B& := fk^k of the field B from fi£ to the 
whole space together with the estimates on the field strength (J5.7|) and on the variation of the 
field direction (|5.8jl . 

Finally, the appropriate vector potential A& is defined as := A + A*, where A* is the 
Poincare gauge of the magnetic field B& — B centered at k. Since B^ = B on Q* k and fll is 
convex, we have A* = on Q£ as well. This completes the proof of Lemma I5TT1 □ 



A. 2 Proof of Lemma 16.11 

Using ()6.1U|) and D n = D 3 , we compute the commutator: 

[D*ZD, D 3 ) = [DjZjkDk, D 3 ) = /;,. D 3 ] + iD^Z^Dk + [D jt D 3 ]Z jk D k 

= (-i)(D j Z jk (Q%D rn -i(d m Q%)+Pa) + ^ 

= (— i) ^DjZjkQk%D m — ^(d m Q'Jl)(Zkj — Zjk)Dk — \dj (Z jkidmQ 1 ^)) + Zjkf3k?jDj — idj(ZjkPkz) 

+ (DmQfa + Pjt) ZjkD k - Dj(dzZ jh )Dkj , 

and 

[h,D 3 ] = id 3 h. 

Now we match the coefficients in 

[D*ZD + h,D 3 ] = 0. 
Matching the highest (second) order terms gives 

d 3 Z = ZQ + Q l Z 

with the matrix Q a b := Q b a3 . 

On the supporting plane we choose Z = I. Using ()6.4j) we have V 7 <5 = 0(£~ 7 L~ 1-7 ), 
< 7 < 4, in a neighborhood of size 0(eL), otherwise Q = 0, hence we obtain the solution Z 
satisfying (JbM4j) . and Z is clearly a symmetric, real matrix. 
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Using the established symmetry of Z, to match the first order terms, we need that 

Z kj (3 j3 = k= 1,2,3 

but clearly j3j 3 = because e 3 = n was the direction of the magnetic field. 

Finally the constant term can be matched by choosing h to be on the supporting plane 
and such that h satisfy 

d 3 h = - 1 1 d J (Z lk (d m QT 3 )). 

Using again that V 7 Q = 0(e~ 7 L~ 1-7 ), < 7 < 4, in a neighborhood of size 0(eL), otherwise 
Q = 0, we obtain that the solution satisfies (jbM5j) . □ 

A. 3 Proof of Lemma 16.21 

Proof of \6.1§) . We use V — cr ■ II, ()6.7|) . the self-adjointness of Dj, ()6.4|) and a Schwarz' 
inequality to obtain 

3 3 
V 2 = (j>^ + fdiv ej + \a ■ o;(e,)]) 2 > ^($>^) 2 ~ " ^ 

3=1 3=1 

Furthermore 

3 2 2 

3=1 i=i 

where <7_i_--D_i_ := o l D\ + a 2 -D2 and {X, Y} = XY + YX denotes the anticommutator. Using 
that aV 3 = — ia 2 and a 2 cr 3 = ia 1 , we obtain from (jfi.lOJ) and /3(e3, ej) = 0, j = 1, 2, that 

{^■^,^3} = ta 2 [D 3 ,D 1 ]-ia 1 [D 3 ,D 2 ] 

2 2 

= £ a a M ab D b + W + I ]T(-1) V a (g| a D 3 + (A.3) 

a, 6=1 a=l 

with M ab := \{.—l) a Q 3 ^ a,b = 1, 2, where 1 := 2 and 2 := 1. We use Schwarz' inequality and 
that Q l - k = OIL" 1 ) to estimate the last terms containing D 3 . In summary, we obtain from 



that 



V 2 > - 



o- ± -D ± Y + ^D 2 + a a M ab D b + D b M ab a a - 0(L~ 2 ) . (A.4) 
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We compute a two dimensional Lichnerowicz formula 



{a ± -D ± ) 2 = D\ + D\ + ia 3 [D u D 2 ] =D\ + D\ + Ba 3 + \a 3 (Q™D m + D m Q 



12 J 



using that (3i2 = B. Let P± := ^(1 ± cr 3 ) and note that P_, P + are orthogonal projections 
commuting with a 3 , and P±cr a = (T a P T for a = 1,2. In particular (cr_i_ -PjJ 2 commutes with 
P±. Therefore 



(<t±-Dj 



P+K-P±) 2 P+ + P-K-P±) 2 P- > P+(^±-P±) 2 P+ 



and 



P+(a ± -D ± ) 2 P + = P + {Dl + Dl + B)P + + \P + a z {Q™D m + D m Q™)P + 

> P + {D\ + D\ + P)P + - \P + {D\ + P 2 + P 3 2 )P + - 0(L~ 2 ) 
by a Schwarz' inequality and ([6.4)1 . We obtain from (|A.4|) that 



p 2 > i 



P + (P 2 + P 2 2 + P)P + + P 2 



0{V 



Finally, we compute 
a a M ab D b + D b M ab a a 



P 



P +P_ 



2P_o a M ab D b P + + 2P + D b M ab a a P- 
+P + a a [M ab , D b ]P_ + P_[P fe , M afe ]<x"P+ . 



a a M afe P b + P b M ab a a 



(A.5) 



P + 



By a Schwarz' inequality and ()6.4|) we obtain 

a a M afe P b + D b M ab a a > ~P+(D\ + P 2 )P+ - O^iT 2 ) 
which, combined with (|A.5|) gives (jOHj) since D 3 = P n . 
Proof of \6.19\i . From the Lichnerowicz formula (j6.13|) and ()6.7|) we have 



D d - | dive,- + \a • u>( ei )J (P; + | dive,- 
D 2 + i[o- ■ u>( ei )P,- + D i0 - • w( ei )] + f , 



u>(e,-)) + <r(*/3) 



where £ is a zeroth order error that satisfies 

\£\ + L\V£\ <(P + b) 

from JO) and IjCT) . 



(A.6) 



(A.7) 
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Lemma A.l. The following comparison estimates hold 



±D 2 - c(P + b) < V 2 < 2D 2 + c(P + b) , (A.8) 

i(D 2 ) 2 - c(P + bf <V 4 < 4(D 2 ) 2 + c(P + bf , (A.9) 
±(D 2 ) 3 - c(P + bf < & < 4(D 2 ) 3 + c(P + 6) 3 . (A.10) 



Proof. The proof of ()A.8|) is straightforward by a Schwarz' inequality and estimating the 
geometric terms as 

( CT -o;( e ,)) 2 <|a;(e,)| 2 <, 2 P (All) 
from ()6.4|) . For ()A.9|) we again use a Schwarz' inequality 



V 4 < 3(D 2 ) 2 + 3 



Djicr ■ a;( ei )) 2 D, + c ■ u{ ej )D 2 jt r ■ u>( ej ) + 3(P + b) 



The first term in the square bracket is bounded by e 2 PD 2 using (jA.llJ) . In the second term 
we commute through and obtain 

(t ■ uj(ej)D 2 cr ■ ufa) < 2P i (<r • u}( ej )) 2 Dj + 2\dj<r ■ uj( ej )\ 2 < 2e 2 PD 2 + e 4 P 2 

from (j6.4j) . Altogether, with a further Schwarz' inequality we obtain the upper bound in (|A.9|) . 
A Similar argument gives a lower bound as well. 

Finally, to bound the sixth power from below we use (|A.8)1 

X> 6 > P 2 (|D 2 — c(P + b))T> 2 . (A.12) 

For the first term we use (|A.6j) and Schwarz 

V 2 B 2 V 2 > i(D 2 ) 3 - ca- ■ uj^D^Dja- ■ u( ej ) - cD j( r ■ cj( ei )D 2 <r • ^{e^Dj - c8T> 2 8 

All the terms cr • oj{ej) can be commuted in the middle and estimated by eP from (|fj.4jl . 
Similarly, Dj's can be commuted through D 2 to estimate 

ePDjtfDj < 2eP(D 2 ) 2 + eP 2 D 2 + eP{P + b) 2 

using ()6.10|) . The resulting terms of the form eP(D 2 ) 2 and eP 2 D 2 can be controlled by 
2e(D 2 ) 3 + ceP 3 . The last term in (|A~T2|) can be dealt with in the same way after estimating 
V 4 by (D 2 ) 2 from (jA~9"jl . For the term ST> 2 8 we use 

£T> 2 £ < 2D J £ 2 D j + 2\\/£\ 2 
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and the bound (jA.7j) . 

Putting all these estimates together and applying several Schwarz' inequality, we easily 
arrive at the lower bound in (jA.10|) . The upper bound can be proven similarly, although we 
will not need it. □ 

After these preparations, now turn to the actual proof of ()6.19|) . 

T> 3 9 2 & > ^V 2 9V 2 9V 2 - 2V 2 \Ve\ 2 V 2 > )S 2 QV 2 0T) 2 - eX&rfV 2 . (A. 13) 

The error can be absorbed into the second term on the left of (j6.19|) . We continue with the 
main term 

V 2 9V 2 eV 2 > -V9V A eV - 6Va(d9)V 2 a(d9)V. (A. 14) 



In the error term we use (jA.8|) and |V#| 2 < ePuj 2 

5Va(d9)V 2 a(d9)V < cSVa(d9)T> 2 a(d9)V + e5P(P + b)Vuo 2 V . (A.15) 

The second term can be absorbed into the third term on the left side of ()6.19|) since 5(P + b) < 
2A and P < A. In the first term in (jA.15|) we can eliminate the cr's using a Schwarz' inequality, 
(dj9j)a j ~D 2 a k (d k 9) < 3(dj9j)~D 2 (dj9j), and the fact that D 2 is identity in the spin space, hence 
it commutes with the Pauli matrices. After that we can again use (|A.8|) and commute V# 
inside 

SVa(de)T> 2 a(d9)V < c5P(9^)D 2 (6» j 0)P 

< cSVidjO) (& + P + b) {djd)V 

< c5V 2 \V9\lV 2 + c8V\ V 2 9\ 2 V + c5eP(P + b)Vuo 2 V 

< ceXV 2 rfV 2 + ce\ 2 Vuj 2 V , (A. 16) 

and these terms can be absorbed into the lower order terms on the left side of (|6.19|) . 
Finally, we continue the estimate of the main term in (jA.14|) : 

-mV^OV > —VQVHV > —dV & d - 5 2 a(d9)V 4 a(d9) . (A. 17) 

In the first term we use flA.10|) after the estimate 5 2 > 5 3 . In the error term we change T> A to 
(D 2 ) 2 using the upper bound in (jA.9|) . then we can eliminate the ex's, use the lower bound in 
(jA.9|) to change (D 2 ) 2 back to T> A and finally commute d8's in the middle similarly to (jA.16|) . 
All the error terms can be absorbed into the lower order terms on the left side of ()6.19|) . We 
omit the details. 
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Putting all these estimates together, we obtain for sufficiently small e 

v 3 e 2 v 3 + xb 2 rfv 2 + x 2 Vu 2 v + x 3 x 2 > c5 3 e(D 2 ) 3 e , 



Since T = D • ZD + h and Z, h satisfy ([6.14)1 . (|6.15|) . in particular Z is close to the identity, 
it is a trivial exercise with commutators to verify that 

(D 2 ) 3 > cT 3 - eP 3 

and the error can be absorbed into A 3 % 2 . This completes the proof of (jbM9j) . □ 
A. 4 Proof of Lemma 16.31 

Since X := 5\~ 1 T and Y := \~ 1 D 2 l commute, we can estimate 



A 3 



< 



(ST) 5 + \ 2 Dl + A 3 - (X + yi/3 + J) 



< c / t 2 e-*e- tx e- iyl/3 dt 



therefore 



x 



(STy + VD i + \> {x - x)ic l iV ' 



-tx. 



e- tYl/ \y,x) 



dy dt . 



(A.18) 



We recall that T and Dj refer to the nonmagnetic counterparts of T and Dj, we also denote 
X := 5A~ 1 T, Y := A _1 D 2 . The operator T is a second order uniformly elliptic operator, its 
heat kernel has a Feynman-Kac representation with an imaginary term in the exponent due 
to a. Following the standard proof of the diamagnetic inequality using Feynman-Kac formula 
(see e.g. (S79j ). we can estimate the oscillatory term trivially by one and we obtain 



-tx. 



x,y) 



< e 



-tx, 



x,y) 



From standard off-diagonal elliptic heat kernel estimates (see, e.g. |LiYj ) we obtain 

c(x - y) 2 ^ 



-tx 



(x,y) < citSX' 1 )- s > ' 2 exp 



tsx- 1 

The operator D n and D n are conjugated by a phase factor 



(A.19) 
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where is a solution of d n q 
exists globally. Therefore 

hence, by spectral theorem, 



a(n). Since the integral curves of n are disjoint, a solution 

Y = e'^Ye 1 ' 1 ' 



e (y,x) 



e (y,x) 



Since Y is a one-dimensional Laplacian, its spectral decomposition is exactly computable 
in appropriate coordinates (£1,^2, £3) chosen as follows. Let £3 be the arclength parameter 
along the integral curves of n with £ 3 = on the supporting plane V. We choose an arbitrary 
Euclidean coordinatization £j_ := (£1, £2) on V. Since the integral curves give a regular foliation 
of R 3 , we can extend the coordinates £1, £2 onto R 3 by setting it constant along each integral 
curve. 

Since the integral curves are straight lines outside of and they are almost parallel 
inside, it is easy to see that 

ie(*)-e(y)i 



Cl < 



\x - y\ 



< C 2 



(A.20) 



for any x, y with universal constants ci, c%. A detailed proof of this obvious fact can be obtained 
along the lines of the proof of Lemma 8.4 of |ES-lVj . therefore we omit it here. Moreover, the 
Jacobian determinant det(D£(x)) of the transformation x 1— > £(x) is universally bounded and 
likewise for the inverse map. 

Define the unitary transform U : L 2 (R 3 , d£) — ► L 2 (R 3 , dx) 

[Ui/>] (a;) = |det(De(x))| 1 / 2 ^(e(x)). 



We then have that 



U*D n U 



-id, 



Hence 



Hence 



x,y) 



det(DH(x))\ 1/2 \det(Dii(y))\ ir2 e 



1/2 -t[(-id is f] 1/3 



(teUG/)) 



det{D£{x))\ 1/2 \ det{D£{y))\V 2 6{Z ± {x),£ ± {y)) / e - tA " 1/s W a/3 e^ 3 ^- £3 ^dp. 



R 



-tyvs, 



x,y) <c\VH-V 2 5(U(x),U(y))- 



(A.2i; 



Now we can estimate the diagonal kernel (|A.18|) . using (|A.19|) . (jA.21|) . (|A.20p . and the 

boundedness of the Jacobian determinants; 



A 3 



(<5T) 3 + X 2 D 2 n + A 3 
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< cX'S- 3 " J" J rV< exp ( - ''^..'^ feiiMlitldl 



:A 3/2 5 _1 . (A.22) 
This completes the proof of (|6.20j) . □ 
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